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It is demonstrated through Monte Carlo simulations that the one component lattice Coulomb gas model 
in two dimensions under certain conditions display features of an anomalous dynamic response. We suggest 
that pinning, which can either be due to the underlying discrete lattice or induced by disorder, is an essential 
ingredient behind this anomalous behavior The results are discussed in relation to other situations where this 
response type appears, in particular the two components neutral Coulomb gas below the Kosterlitz-Thouless 
transition, as well as in relation to other findings from theory, simulations, and experiments on superconductors. 
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I. INTRODUCTION 



Two-dimensional (2D) vortex physics is strongly reflected 
in the properties of systems related to 2D superfluids like 
Josepjhson junction arrays,El'a superconductin|aiilms, and '*He 
filmsjj as well as high-Tc superconductors .Hcl The vortices 
are, e.g., respotmhle for the well-known Kosterlitz-Thouless 
(KT) transitionElu'tl and there is a fairly good understanding 
of the thermodynamic properties related to 2D vortex physics. 
However, the dynamical aspects, which can be probed by the 
complex imcedance and the flux noise measurements in super- 
conductorsJl BEl and by the torsional oscillator period shift in 
^He-fiIms,Q are much less well understood. The first attempt 
to describe the dynamic vortex response phenomenologically 
was given by Ambegaokar-Halperin-Nelson-Siggia (AHNS) 
in Ref. ^ However, it was later discovered that the Minnhagen 
phenomenology (MP) in Ref. |l] based on essentially the same 
ingredieute as AHNS, described experiments and simulations 
betterLrt£l A particularly interesting aspect in the. MP is that 
it reflects an anomaly in the 2D vortex responsejlj which, for 
example, in a superconductor takes the form of a logarithmic 
divergence of the complex conductivity: cr(cj) cx — Incj for 
small frequency lo. There still remain open questions on what 
ingredients can cause this anomaly and give rise to the anoma- 
lous response form. 

There has been a number of attempts to obtain the MP form 
with other phenomenologicalj^ipfpaches as well as somewhat 
formally more rigorous ones.EJTH3 One possibility is that it is 
an intrinsic property of the 2D vortex system which is linked 
to the long range logarithmic interaction, as presumed by the 
original motivation for the MP response form.H Another pos- 
sibility is that it is caused by the coupling between the perpen- 
dicular currents described by the vortices and the Ipjpgitudinal 
currents (often referred to as the spin-wave part).0'E3 In the 
present paper we suggest that there is yet another possibility; 
the MP type response can also arise from pinning influencing 
the vortex system. 

The approach to this problem, which we subscribe to here, 
is to systematically investigate various models by aid of sim- 
ulations. Vortices can be described as 2D Coulomb gas parti- 
cles, in the sense that a vortex with vorticity ±1 corresponds 
to a Coulomb gas particle with charge ±1 (See, e.g., Ref. |l] 
for details on the mapping between the two models). Since 



the dynamics of the vortices can to a |g|E)od approximation be 
described by the Langevin dynamics,El B we can view the 2D 
Coulomb gas model with Langevin dynamics as the dynamic 
model for the "pure" 2D vortex response which includes nei- 
ther any coupling to a spin-wave part nor any pinning. For a 
superconducting film in the absence of an external perpendic- 
ular magnetic field (which corresponds to the absence of a net 
rotation in case of superfluid ^He film), the numbers of vor- 
tices and antivortices (with vorticities 1 and —1, respectively) 
are always equal, and the system is described byj-the charge 
neutral two components 2D Coulomb gas modeljil on which 
the AHNS, as well as the MP, are based. Although it was 
clearly shown from simulation of this neutral 2D Coulomb 
gas model that the low-temperature phase is well described 
by the MP respopse form with the logarithmic divergence of 
the conductivityjlj there still remains the possibility that the 
coupling to the spin-wave part or pinning could also influence 
and perhaps lead to the anomalous response in certain cases. 

In Ref. [ri| it has been found in experiments that the 2D tri- 
angular Josephson array in a perpendicular magnetic field is 
also well described by the MP response form, which has again 
been confirmed in simulations ofJ±ej:eIated frustrated 2D XY 
model with relaxation dynamicsollj Since in these cases pin- 
ning could with large probability be excluded as the cause of 
the effect, the observed anomalous behavior can in these cases 
be attributed to either the logarithmic vortex interaction, or to 
the interplay between the vortex and the jsain-wave interac- 
tions, or to a combination of both.EilEJOilj From the analy- 
sis of the experimental results in Ref. ^ for thin MoGe super- 
conducting film in a perpendicular magnetic field.|_Lt has also 
been pointed out that the data show MP behavior.Ej Since in 
Ref. ^ it was argued that for this sample the vortex pinning 
due to disorder was important, this raises the question of how 
pinning influences the dynamic behavior Could the MP re- 
sponse form also arise from an interplay between the vortex 
system and pinning? In the present work we conclude, from 
simulations of the lattice Coulomb gas model, that the vortex 
system, in the absence of couplings to the spin-wave part but 
in the presence of pinning, indeed in certain parameter regions 
gives rise to MP features reflecting an anomalous diffusion. 

The paper is organized as follows: In Sec. ^ we describe the 
Coulomb gas model and the correlation functions we study in 
simulations. Section |ll| gives a brief description of the fea- 
tures of the MP response form and make comparisons with 
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the standard Drude response form. The resuhs from the simu- 
lations are given in Sec. |^ and in Sec. we try to give some 
perspective from comparisons with other results. 



II. THE LATTICE COULOMB GAS MODEL 

In this work, we use the 2D Coulomb gas model on an L x L 
triangular lattic&qwhose Hamiltonian in the absence of disor- 
der is written agHJ 



(1) 



where n, is the integer charge on the ith site and the frustration 
/ controls the total number of charges Nc = J2i = 
{N = L^). The interaction Vij between charges at positions 
r i and Vj for a triangular lattice is given by 
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N ^ 6 — 2 cos k • ai — 2 cos k • ao — 2 cos k • a-? ' 
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where ai = a, a2 = b, and aa = a — b with two primitive 
translational vectors a and b for 2D triangular lattice. (We 
choose a = X and b = x/2 + \/3y/2 with lattice spacing set 
to unity.) 

We also study the effects of quenched disorder and consider 
the Hamiltonian 



(3) 



where the pinning potential Ui at i has value C/p(> 0) if 
the disorder is located at i and J7i = otherwise. The 
randomly distributed disorder realization is parameterized by 
the pinning strength and the ratio r.p between the num- 
ber of pinning sites and the total number of charges N^'- 
rp = Np/Nc- This corresponds to point disorder since we can 

express it as U{y) = X^jl^i UP5{r — Vj) in continuum limit, 
where rj denotes the positions of disordered sites. 

We use Monte Carlo (MC) dynamics which can be imple- 
mented as follows: Allow two different MC tries where one 
is the charge hopping to one of its nearest neighbors, i.e., 
(1, 0) (0, 1), and the other is the generation of a charge- 
anticharge dipole on one bond, i.e., (0,0) (1,-1). As 
was already observed in Ref. ^ for parameters used in this 
work the acceptance ratio for generation of charge dipole is 
found to be extremely small even at temperatures much higher 
than the critical temperatures (referring to the depinning and 
melting of the solid phase). We are here studying a temper- 
ature region where no charge dipoles are in practice created. 
Consequently, we only include the charge hopping in our MC 
try. This means that the model we are studying only contain 
charges of one sign. Our algorithm is thus: At each step, we 
first pick up one charge at random and try to move it to one of 



its six nearest neighboring sites which is also randomly cho- 
sen. This MC try is accepted or not according to the standard 
Metropolis algorithm. One can make the MC simulation more 
efficient by using a nonlocal update (e.g., the randomly chosen 
charge can be allowed to hop to a distant site). On the other 
hand, as far as the dynamic behavior is concerned, the local 
update rule like the one used in the present work is believed 
to better imitate the relaxational dynamics of the correspond- 
ing continuum model (which corresponds to the Coulomb gas 
model with Langevin dynamics).Ej 

For the characterizations of the thermodynamic properties, 
we measure (i) the dielectric constant l/e(0), which corre- 
sponds to the helicity modulus in the vortex systems and de- 
tects the superconducting to normal transition, (ii) the orien- 
tational order parameter (f)^, which probes the melting of the 
solid-like structure, and (iii) the specific heat Cy (see Ref. e3): 
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where (• • •) is the thermal average, = n^e"*'^ '"' is the 
Fourier transformation of the charge distribution, Oi is the an- 
gle between a reference direction, say x, and the line con- 
necting the ith particle and its closest neighbor. Since the 
smallest wavevector is limited by the system size, we choose 
k = A/L, B/L, and -(A + B)/L with reciprocal lat- 
tice vectors A = 47ry/\/3 and B = 27r(x — y/\/3), and 
take the average over these three smallest k's (with magni- 
tude k = y/l6TT^/3N) to obtain the static dielectric function 
l/e(0) in Eq. (§. 

The dynamic behaviors are described by the dynamic di- 
electric function 1 /e(ti') (see Refs. n2 and 13 for details): 
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G(t) = lim — 
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where (?T.k(i)?T._k(0)) is the charge correlation function, and 
l/e(0) is the static dielectric constant given above. Just as 
for l/e(0) in Eq. (Q), G{t) is obtained as the average over the 
three smallest wavevectors. We measure time t in units of the 
MC step so that one time unit corresponds to Nc MC tries to 
move the randomly chosen particle to one of its nearest neigh- 
bors. 



III. DYNAMIC RESPONSE 

The dynamic response function we focus on is the complex 
dielectric function l/e{uj) given by Eq. (R). The anomalous 
vpitex diffusion is reflected in the fact Gjl) oc l/t for large 
tBB From this one directly infers that to leading order in small 
uj (Ref. |l6|) 
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where 1/e is a constant and the first equality in Eq. (10) fol 
lows because l/e(0) is always a real quantity. The com 
plex conductivity for a superconductor a{Lo) corresponds t| 
— \/ioje{Ld) which means that Yiea{uo) cx — hiu) for small lo 
This logarithmic divergence is another characteristics of the 
anomalous vortex diffusion. 
The MP form is given by 
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and the corresponding explicit form of G{t), which we de- 
note by Gmp, hence incorporates the features associated 
with G{t) oc l/tJl3 The MP form is associated with a sin- 
gle characteristic frequency ojq. In the 2D neutral Coulomb 
gas all vortices are boiipjd in vortex-antivortex pairs in the 
low-temperature phaseE'Q This implies an infinite correlation 
length which for a superconductor means a vanishing resis- 
tance Ra Since the logarithmic divergence of the conductiv- 
ity means that i? = 0, a possible scenario is that the MP re- 
sponse form describes the response of the vortex-antivortex 
bound pairs and hence should apply to the low-temperature 
phase. It has been verified, from simulations of the 2D neu- 
tral Coulomb gas with Langevin dynamics, that the MP form 
indeed gives a very good description in this case.EJ Above the 
KT transition there are both free vortices and bound vortex- 
antivortex pairs present and hence there is a finite correlation 
length. This means that G{t) has an exponential decay of the 
form e^*/'", where t is a relaxation time and i? oc I/t. Thus 
a finite resistance is associated with an exponential decay of 
G{t). If there were only free vortices then G{t) would be 
dominated by the exponential factor and the response should 
to good approximation be described by the standard Drude 
form Gd oc e^*/'^, i.e.. 
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More generally, when free vortices are present the response 
form Gmpd = GMp(O^D(i) has from simulatinns-ieen 
shown to give a good parameterization of the data.Ljllj€Zl 

One difference between the MP response and the Drude 
response, which is easy to gauge in experiments and simu- 
lations, is the peak ratio: The peak ratio (PR) is defined as 
the ratio between the real and imaginary part of the com- 
plex dielectric function at the frequency where |Im[l/e(a;)]| 
has its maximum, i.e., at the dissipation peak (see Fig. 6). 
For the standard Drude response this ratio is 1 whereas it is 



2/tt K, 0.64 for the MP response. In general when the re- 
sistance is finite the peak ratio should be between these two 
values, 2/7r < PR < 1. 

Our strategy is to calculate G{t) as described above and 
then to gauge the response by studying the peak ratio, G{t) 
for large t, and to what extent the response is parameterized 
by the MP form. 



IV. SIMULATION RESULTS 



A. Lattice Coulomb gas without disorder 



In this section, we first present the results of the MC simula- 
tions of the lattice Coulomb gas model given by the Hamilto- 
nian in Eq. ([I]) for the case without random pinning sites, i.e., 
Ui = Q \n Eq. (|]). It is known that the system can undergo 
two phase transitions in this case:EJ One is the depinning of 
the solid at Tc, corresponding to the superconducting to nor- 
mal transition, and the other is the solid to liquid transition at 
Tm (Tc < Tm)- These two transitions can be identified by 
aid of the quantities l/e(0) and 4)q, respectively. Since the 
continuum system corresponds to an infinitesimally small lat- 
tice spacing and the frustration / is proportional to the area 
of the elementary plaquette, the continuum limit of the lattice 
Coulomb gas is given by / 0, but with the number of par- 
ticles Nc still finite: Nc = A^/=const. The observation in 
Ref. ^ that Tc becomes smaller as / is decreased indicates 
that Tc ^ in the continuum system so that the system only 
has one transition at a nonzero temperature Tm- 

We start by presenting some static results for / = 1/16 
on a 64 X 64 triangular grid. In the simulations, we start 
from high enough temperatures and anneal the system by de- 
creasing the temperature slowly, to avoid being captured by 
metastable states at low temperatures. However, the local MC 
scheme used in this work (a particle is only allowed to hop to 
its nearest neighbor sites) makes it difficult to achieve equi- 
libration near and below the transition. Thus very many MC 
updates are required in this region (our longest runs consists 
of 2 X 10'5 MC steps). As T is decreased l/e(0) abruptly 
changes to unity near T « 0.01 (see Fig. |l]). This abrupt 
change is associated with the normal (high T) to supercon- 
ducting (low T) transition in the vortex system since l/e(0) 
corresponds to the superfluid densityE In the Coulomb gas 
model it corresponds to a transition where the Coulomb gas 
systems becomes pinned to the underlying grid. The six- 
fold orientational order parameter 06 has the value unity for a 
perfect triangular Abrikosov vortex lattice. Consequently the 
abrupt drop of 06 seen in Fig. |l] indicates that the vortex solid 
melts into a vortex liquid near T w 0.01. The specific heat G„ 
also shows a sharp maximum near this temperature T k, 0.01. 
Consequently the data in Fig. |l| for the lattice Coulomb gas 
model with / — 1/16 are consistent with a single phase tran- 
sition corresponding to a transition where the superconduc- 
tor becomes normal apd the Abrikosov vortex structure melts, 
i.e., Tc ~ T„, » O.OlEil 
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We next turn to the dynamic results for / = 1/16. These 
are shown in Fig. ^ The characteristic frequency ujq in 
Fig. |2|(a), determined from the peak position of |Im[l /e(aj)] |, 
shows a dip-like feature near T = 0.01, reflecting a critical 
slowing down close to the phase transition (compare the static 
results in Fig. |l|). In general, a system composed of particles 
is expected to be well described by a Drude approximation at 
sufficiently high temperatures because the interaction between 
the particles becomes negligible in comparison with the strong 
thermal fluctuations. The dynamic response function in such 
a temperature regime should be well described by the Drude 
response form Eqs. and (|o|) and should consequently be 
well characterized by the peak ratio (PR) equal to one. As 
is shown in Fig. |l](b), the dynamic dielectric function indeed 
approaches the Drude PR value one at higher temperatures. 
However, the PR crosses over to a value close to PR=0.64 
which characterizes the anomalous behavior (sea the previous 
section) as we approach T « 0.01 from above.a 

The question is then what causes this crossover behavior. 
Is it associated with the depinning transition at Tc or with the 
melting transition at Tm? In order to address this question one 
needs a temperature separation between the two transitions. In 
Ref. ^ it was shown that as / is decreased, the two transi- 
tions present become well separated. With this in mind we, as 
a next step, repeat the simulations for a much smaller value of 
/• 

Figure |] shows the static results l/e(0), 06, and Cy for the 
system with / = 1/64 on a 96 x 96 lattice. As seen from 
the figure, the superconducting to normal transition is now at 
Tc « 0.0035 and is well separated from the vortex solid melt- 
ing transition at T„i w 0.007, as was found in Ref. As also 
seen in Fig ||, the existence of two transitions are also reflected 
by two peaks in C„: one near Tc and another broader one near 
Tm- The intermediate phase existing between Tc and Tm is 
characterized by a vanishing l/e(0) together with a nonzero 
(fiQ and is interpreted as a floating solid phaserJ The vortex 
solid structure is depinned from the underlying discrete grid 
and floats around, which leads to dissipation of energy. 

Next we analyze the dynamic response function for / = 
1/64 is terms of the characteristic frequency ujq and the peak 
ratio PR. In Fig. ^ it is shown that luq displays a dip-like struc- 
ture-Jjeflecting a critical slowing down near Tc rather than at 
Tm.Ej Similarly, the correspondence of the crossover behav- 
ior for / = 1/16 in PR from Drude to anomalous response 
form shown in Fi*. ^b), is in Fig. |(b) found to be near Tc 
and not near T,„.cJ From this we conclude that the crossover 
behavior in the dynamic response is associated with Tc and 
hence with the pinning of the Coulomb gas to the underlying 
lattice. Thus for the pure one component lattice Coulomb gas 
with MC dynamics we infer that there is a crossover behavior 
from a normal Drude like response to an anomalous response. 
An essential ingredient for the appearance of the anomalous 
response is associated with the pinning between the vortex 
system and the underlying lattice. 



B. Lattice Coulomb gas with disorder 

In the previous section we inferred that the pinning, caused 
by the underlying lattice, leads to a crossover from a Drude 
to an anomalous dynamic response. In the present section we 
study the role of pinning further by introducing extrinsic vor- 
tex pinning caused by disorder into the model. The general 
idea is that the anomaly in the dynamics is caused by a slug- 
gish motion of the vortex system and that pinning hindering 
the vortex motion can cause such a behavior. To this end we 
introduce point impurities into the lattice Coulomb gas model 
as described by the Hamiltonian in Eq. (|^). For a supercon- 
ductor this corresponds to magnetic point impurities. In this 
section, we perform MC dynamic simulations for 10 different 
disorder realizations in order to get disorder averaged quanti- 
ties. 

In Fig. H the PR for the system with / = 1 /64 and L = 96 
atT ~ 0.01 is displayed as a function of the pinning strength 
JJP (the number of pinning sites Np is fixed to 14, which ap- 
proximately corresponds to 10% of the total number of vor- 
tices Nc, i.e., Np w O.lA^'c). As seen in Fig. || the dynamic 
response for the pure case {U^ ~ 0) has PR close to one at 
this higher T [compare Fig. 0(b)]. This means that the pin- 
ning to the underlying lattice at this T is too weak to influ- 
ence the behavior so that the vortices should obey a Drude re- 
sponse. However, if we introduce disorder then Fig. |] shows 
a decreasing PR with increasing pinning strength again sug- 
gesting a crossover towards an anomalous dynamic response 
[compare Fig. ^b)]. This is consistent with the idea that when 
the pinning strength becomes stronger, the motion of the vor- 
tices becomes more sluggish, resulting in a decrease of the PR 
and a crossover to an anomalous response. 

This scenario implies that if we introduce enough pinning 
then the crossover to the anomalous response should become 
complete. To achieve this we investigate a stronger pinning 
case where Np = Nc and the pinning strength = 0.1. 
Figure ^ shows the obtained dielectric function l/e(w) for 
this case [see Eqs. (Q) and (^] at T = 0.02 which is much 
higher than Tc as well as T„i for the pure case (Tc ~ 0.0035 
and Tm ~ 0.007, respectively). As illustrated in Fig. ||, the 
PR is indeed close to the anomalous value PR= 2/7r in this 
case. In the MP phenomenology the anomalous dynamics is 
linked to a 1/i decay of the correlation function G{t) [com- 
pare Eqs. (p|)-(p^]. The time-correlation function G{t) (cor- 
responding to 1 /e{uj) in Fig. ^ is shown in Fig. ^ and indeed 
has a 1/t-tail as is manifested by the horizontal plateau of 
tG{t) for larger t. This suggests that there exists a link be- 
tween a 1/t-tail in G{t) and the PR= 2/tt. Such a link is 
incorporated into the MP response form Eqs. (^l]) and (|l^). 
The inset in Fig. ^ shows what happens for the same case at 
the higher T = 0.03. In this case G{t) suggests the pres- 
ence of an exponential decay. This is very similar to previous 
simulations where Gmpd = GMp{t)GD{t) has been shown 
to give a good parameterization [see Sec. |l| below Eq. ( p4| ) 
and Refs. |l3[ Q and ^7|] . Figure ^a) shows that a reasonable 
fit of l/e{ui) to the MP form is obtained in the anomalous 
case T = 0.02. This suggests that the basic link between 
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the G{t) (x 1/t and the PR=2/7r is reasonably well captured 
by the MP form. Figure ^b) shows a fit to the exponen- 
tially decaying case at T = 0.03 using the parameterization 
Gmpd = GMp{t)GD{t). Again a reasonable fit is obtained 
with PR« 0.9. This suggests that, although it is close to the 
Drude form, there is a small crossover towards the anomalous 
response (see Fig. ||). 

In this section we have thus shown that for the one compo- 
nent lattice Coulomb gas, at a T so high that the pinning to the 
underlying lattice plays no role, it is still possible to obtain a 
crossover to the anomalous response by introducing pinning 
through randomly distributed pinning sites. 



V. COMPARISON AND DISCUSSIONS 

We have found from simulations that the 2D one component 
lattice Coulomb gas with Monte Carlo dynamics under cer- 
tain circumstances displays an anomalous dynamic response. 
An essential ingredient for the appearance of this anomalous 
response is found to be pinning; either pinning due to the un- 
derlying lattice or external pinning induced by disorder The 
anomalous response vanishes gradually as the pinning effects 
become weaker and from this we conclude that there is no 
anomaly in the absence of pinning. 

We gauge the degree of anomaly in the response by the be- 
havior of the correlation function G{t) for large t, the peak 
ratio, and to what extent the response can be parameterized by 
the MP form. We find that the full anomaly is characterized 
by G{t) cx 1/t for large t, a peak ratio consistent with 2/7r 
and note that these two features are also incorporated into the 
phenomenological MP form. For the one component lattice 
Coulomb gas this means that an anomalous response in the 
absence of externally introduced pinning only appears very 
close to the depinning of the vortex system from the underly- 
ing lattice structure. From this we further infer that the pure 
continuum one component Coulomb gas has no anomalous 
response and that in this case external pinning is required to 
change the response from normal to anomalous. 

It is interesting to compare this with the 2D neutral (two 
components) continuum Coulomb gas with Langevin dynam- 
ics which is the generic model for the vortex physics of a 2D 
superconductor(superfluid) in the absence of an external per- 
pendicular magnetic field (net rotation).El This model under- 
goes a KT transition and below this transition the dynamic 
response has from simulations been shown to be anomalous 
characterized by G{t) oc 1/t for large t, a peak ratio consis-. 
tent with 2/7r, and to be well described by the MP form.EJ 
Since the same response form is found for the one component 
Coulomb gas in the presence of external pinning, as for the 
2D neutral (two component) Coulomb gas below the KT tran- 
sition, one might speculate about a common physical feature. 
We here suggest that the common physical feature is pinning; 
in the one case externally introduced and in the other intrinsic 
in the sense that negative vortices act as (moving) pinning cen- 
ters for positive ones (and vice versa). In both cases the pin- 
ning mechanism suppresses the free vortex diffusion. When 



this suppression is complete the response becomes anomalous 
as for the one component Coulomb gas with enough exter- 
nal pinning and for the neutral Coulomb gas below the KT- 
transition where all positive and negative vortices are bound 
together in pairs. 

Thus our conclusion based on simulations is that the 2D one 
component continuum Coulomb gas cannot have an anoma- 
lous response by itself. This is in contrast to the theoreti- 
cal considerations in Ref. [l^, based on a Mori approximation 
scheme, which suggest that the one component model in cer- 
tain parameter regions could have an anomalous response. We 
have not been able to find any support for this suggestion in 
our simulations. _ 

The MP formlll was originally motivated for the 2D neu- 
tral (two components) Coulomb gas as an extension of the 
AHNS phenomenologyQ and was assumed to describe the re- 
sponse of bound vortex pairs. Other recent improvements of 
the AWJ|S phenomenology also gives a peak ratio close to 
2 / TT.CJ Ell However, the MP form is to our knowledge the only 
form that incorporates both the large t behavior G{t) cx 1/t 
and the peak ratio 2/7r, both of which seems to be charac- 
teristics of the anomalous response. In the one component 
Coulomb gas case with external pinning there is at present re- 
ally no motivation for the MP form; it is just a simple form 
which simultaneously incorporates two features-which seem 
to be characteristics of the anomalous response.E3 

Our basic conclusion is that the anomalous response form 
can indeed arise from an interplay between pinning and a vor- 
tex system associated with an external perpendicular magnetic 
field. In this context one may note that ^^"^ peak ratio was 
found for the experimental data in Ref. ^llJ This was a MoGe 
superconducting film in a perpendiculaonagnetic field where 
pinning was expected to be important.Ej We suggest that this 
might be an experimental example of our present findings 
from simulations. 
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In connection with Fig. M for the disordered case the possibility of 
a link between PR=2/7r and G{t) ~ 1/f is discussed. This link 
also seems to be borne out in the present case without disorder 
although the convergence of the data is not as impressive as in 
Fig.|. 

The fact that loo does not reach zero in this case is due to the finite 
lattice size. We believe that this is also the reason why the PR does 
not decrease all the way to 2/7r. 

The fact that disorder may cause a slow decay of the correlations 
is reminiscent of glassy states. However, the special feature of the 
anomalous response discussed here is that the correlations decay 
precisely as 1/t. 

The cause of the increase of PR for decreasing T at the lowest 
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FIG. 1. Static dielectric function l/e(<^ = 0), the orientational 
order parameter <^6, and the specific heat Cv vs temperature T for 
the 2D lattice Coulomb gas model on a 64 x 64 triangular grid with 
the frustration / = 1/16. The depinning transition at Tc, corre- 
sponding to the superconducting to normal transition, is signaled by 
an abrupt change of l/e(0). The solid to liquid transition at Tm 
probed by occurs at approximately the same temperature, i.e., 
Tc ~ Tra ~ 0.01. also shows a sharp maximum near the same 
temperature. 
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FIG. 3. Static dielectric function l/e(t<-' = 0), the orientational 
order parameter 06, and the specific heat Cv vs temperature T for 
the system with / = 1/64 and L = 96. The depinning transition is 
at Tc ~ 0.0035 and melting transition is at Tm ~ 0.007. Also note 
that Cv displays two peaks: one near Tc and the other one near Tm.. 
The intermediate phase existing between Tc and Tm has a nonzero 
orientational order and corresponds to a floating solid phase. 
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FIG. 2. Characteristic features of the dynamic dielectric func- 
tion for the system with the frustration / = 1/16 and size L — 64. 
(a) The characteristic frequency scale ljo as a function of T obtained 
from the peak position of |Im[l /e((j;)] | : ojo decreases near T ~ 0.01, 
reflecting a critical slowing down, in accordance with Fig. |lj (b) The 
peak ratio PR defined by the ratio between real and imaginary part 
of l/e(ci;): The response function is described by the simple Drude 
response form with PR= 1 at high enough temperatures and crosses 
over to the anomalous behavior with PR= 2/7r as T is decreased.Lil 
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FIG. 4. Characterization of the dynamic dielectric function for 
the system with the frustration / = 1/64 and size L — 96. (a) 
The characteristic frequency ujo vs temperature T shows a dip-like 
structure near T ^ Tc ~ 0.0035 (see Fig. ^ for the static results), 
(b) The peak ratio PR vs T reflects a crossover behavior from Druda 
response form in high temperatures to the anomalous one near TcBj 
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FIG. 5. The peak ratio PR vs the pinning strength for the 
system with / = 1/64 and i = 96 at T = 0.01. The number of 
pinning sites was fixed to 14 (which is about 10% of the total number 
of vortices and the results correspond to an average over 10 random 
pinning realizations). The PR starts from unity for the pure case 
(U^ — 0), which is the value for the Drude response, and decreases 
as is increased, reflecting that the vortex motion becomes more 
sluggish. 
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FIG. 7. The time-correlation function G(t) for the same case as 
in Fig. ^plotted as In tG{t) vs t. The horizontal plateau for larger t 
shows that the long-time behavior of G{t) has al/t decay. The inset 
shows the same case at T = 0.03 and at this higher temperature G(t) 
appears to have an exponential decay. 
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FIG. 6. The dynamic dielectric function obtained for the system 
with / — 1/64 and L = 96 in the presence of extrinsic pinning 
with = 0.1 and Afp = Nc foi T = 0.02 (10 random pinning 
configurations were used). Filled circles and open circles correspond 
to Re[l/e(cij)] and |Im[l/e(tj)]|. The PR (the ratio between the real 
and imaginary part at the maximum of the imaginary part) is found 
to be close to the MP value 2/tt. 
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FIG. 8. Fits of the dynamic dielectric function to the MP form, (a 
l/e{Lo) corresponding to the anomalous G{t) for T = 0.02 in Fig 
is fitted to Eqs. (jlT]) and ([l2|). Filled circles, open circles and full 
drawn curves correspond to Re[l/e([i;)], |Im[l/e([i;)]|, and the MP 
form, respectively. A reasonably good fit is obtained, (b) l/e{uj), 
corresponding to the exponentially decaying G{t) for T — 0.03 in 
the inset in Fig. ^, is fitted to the parameterization obtained from 
G'MPD(i) = GMp{t)G-D{t) (see the text). Filled circles, open cir- 
cles and full drawn curves correspond to Re[l/e(ti;)], |Iml/e(ci;)|, 
and the l/e(ijj) from Gmpd , respectively. A reasonably good fit is 
found for PR~ 0.9 which is slightly smaller than PR= 1 correspond- 
ing to the pure Drude value. 
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